Gauge symmetry and Slavnov- Taylor identities for randomly stirred fluids 
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The path integral for randomly forced incompressible fluids is shown to have an underlying Becchi- 
Rouet-Stora (BRS) symmetry as a consequence of Galilean invariance. This symmetry must be 
respected to have a consistent generating functional, free from both an overall infinite factor and 
spurious relations amongst correlation functions. We present a procedure for respecting this BRS 
symmetry, akin to gauge fixing in quantum field theory. Relations are derived between correlation 
functions of this gauge fixed, BRS symmetric theory, analogous to the Slavnov- Taylor identities of 
quantum field theory. 
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To be consistent with Newtonian physics, the descrip- 
tion of fluid dynamics and turbulence must be the same 
in all incrtial reference frames. For this reason, the role 
that Galilean invariance plays in the dynamics of a fluid 
system governed by the Navier-Stokes equation has been 
a major focus of study for decades jl|, 01 • The Galilean 
invariance of the Navier-Stokes equation plays an impor- 
tant part in the physical aspects of turbulence [3|, in 
practical aspects of turbulence modelling [4t\ and even in 
fluid animation and computer graphics applications [5|. 
Among the symmetries of the Navier-Stokes equation, 
Galilean invariance is by far the most studied in theo- 
retical approaches to turbulence [l|, S H 0i B 9 ■ The 
physical constraint of Galilean invariance must be taken 
into account in the modelling of subgrid-scale stresses in 
large eddy simulations Q , in models of multi-component 
turbulence [8| , in nonperturbative renormalization group 
analyses of fully developed turbulence 10], and in prob- 
ability density transport equations jTl[. Much of the 
interest has dealt with the infinite number of exact re- 
lations between different correlation functions implied 
by Galilean invariance, which are akin to the Ward- 
Takahashi identities of quantum field theory [l^. The 
most well known of these exact identities relates the ver- 
tex and response functions [21 . Based on this relation, 
various inferences have been made about the nonrenor- 
malization of the advective or inertial term in the Navier- 
Stokes equation [J i, 0, i, 0, [H, [ll ■ 

In this Letter we observe that the Galilean invariant 
generating functional for stirred fluids 2, 15, 16], from 
which all the Ward identities are derived, possesses an 
overall infinite factor, which must first be extracted be- 
fore this functional is well-defined. We draw an anal- 
ogy with abelian gauge theory in quantum field theory, 
whereby the Galilean invariance is associated with gauge 
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invariance, and so the infinity arises from integrating the 
fiuid velocity over gauge equivalent copies (i.e., over an 
infinite number of independent inertial reference frames) 
of the same theory. To remove the infinite factor in the 
Navier-Stokes generating functional, the Fadeev-Popov 
procedure of quantum field theory 1?!] is used to "gauge 
fix" the theory, that is, to fix the theory to a single in- 
ertial reference frame. The gauge fixing results in the 
theory no longer appearing Galilean invariant. How- 
ever, we demonstrate that this gauge-fixed theory pos- 
sesses an underlying Becchi-Rouet-Stora (BRS) symme- 
try [i2, Hj lly] ) which when invoked, restores the Galilean 
invariance. 

This work provides a new insight into the structure of 
the Navier-Stokes equation. However, it is not simply 
a new way to express the same theory. We show that 
the infinite factor associated with the standard dynamic 
functional has associated with it, much more dire conse- 
quences, in that there exist spurious relations amongst 
correlation functions. Explicit examples of such spuri- 
ous relations are given. We then demonstrate that the 
gauge-fixing BRS procedure implemented here not only 
removes the infinite factor, but also eliminates these spu- 
rious relations, thereby rendering a well defined generat- 
ing functional for this theory. The BRS invariant func- 
tional derived in this Letter is the consistent and cor- 
rect functional for the Navier-Stokes theory, that must 
be used to obtain reliable results. 

The standard dynamic generating functional for a 
stirred incompressible fluid has been extensively studied 
for many years @, 0, M, [l3, 14 , [la [la] . It is based on 
the path integral approach ISJ . [2ll . [22 . [23| to classical 
statistical dynamics 24] and is given by 

Z= f[DV][D(T]exp{~S[V,(T]+ f dkdcj{J-V+'E-(T)}, 

where the action 

S'JV.cr] = ^ / dkdui ai{-k, ~uj)Dij{~'k)aj{k,uj) - 



i / dkdu! CTa(— k, —uj) [( — ZCJ + l'k'^')Va{'k,Uj) — 



Mat3j0i.) / didujiVp{k-luj-uJi)V^{i,Ui) 



(2) 



The fluid velocity is V, the conjugate field is or; J and 
S are k and cu dependent sources. Dij{—'k) is the only 
non- vanishing cumulant of the random force, Mai3.y{'k) 
(2i)-i{fc/3Pa7(k) + /c^Pa/3(k)} and Pa^(k) = 6c,p- 






and V is the bare fluid viscosity. Throughout this Letter 
we will work in wavevector-frequency space. 

Consider a second primed reference frame moving with 
respect to the unprimed frame at constant velocity c = 
en where n is a unit vector in the direction of motion. 
The relations between wavevector, and the frequency of 
an event in the two frames and the transformations of 
the velocity and conjugate fields are given by: 

k' = k, w' = w - c • k, (3) 

V(k,w) = V'(k,tj-c-k) + c<53(k)<5(w), (4) 
o-(k,w) = cr'(k,w-c-k). (5) 

The action S and functional measure are invariant under 
this Galilean transformation. The Fadeev-Popov proce- 
dure of quantum field theory [121 Il7l | is now applied to 
show that for J = S = 0, the functional in Eq. ^ has 
an overall infinite factor. For this, expressing the velocity 
zero mode as V(0, 0) = Vo(53(O)^(0) ^, first insert the 
factor of unity 1 — j dh (5'^(Vo — b), into Eq. ([1]), to give 

Z= f[DV][D(T] f dhS^(Vn-h)cxp{-S[V,(T]}. (6) 

The observation is now made that the integrand in Eq. 
([S]) is independent of b. To see this, consider changing 
from b -^■ b + a. Then perform a Galilean transforma- 
tion with boost velocity c = a in Eqs. ([3]) - (O. The 
action and measure are invariant to this transformation 
and (5'^(Vo — (b + a)) — » (5'^(Vo — b), thus restoring the 
integrand to its original form. Due to the independence 
of the integrand on the value of b, the dh integral can 
completely factor out as an overall infinity as 



Z = 



Jdh\ J[DV][DcT]S^{Vo^h)exp{^S[V,(T]}. 



(7) 

The first term in brackets on the right-hand side is the in- 
finite factor we originally sought to isolate. At this stage 
one could simply eliminate this integral since it is just 
an overall, albeit infinite, normalization factor. In this 
case different choices of b, correspond to expressing the 
theory in different inertial reference frames. On the other 
hand, since the integrand in Eq. ([6]) is independent of b, 
we can also insert any function G{h) into the integrand, 
for example to render the b integration finite. One con- 
venient choice is, G(b) — exp(— ^b • b), where ^ > 0. 
Inserting this into Eq. ^ and now performing the dh 
integration, leaves the gauge-fixed (GF) functional 



Z'gf = l[DV][Da]exp{-S[V,a] - ^Vq ■ Vq}. 



In our analogy to gauge theories, all these different pos- 
sibihties for G{h) above would be called gauge choices. 

Problems also arise if one attempts to compute any 
correlation function using Eq. ([T]). For example, consider 
(VQ(k, a;)V^(k', oj')). Using the Galilean invariance of the 
action, this can be written as 

( Vo.{k,Lj)V0{-k',Lo'))= f[DV][D(T]cxp{-S[V,cT]} 



(9) 



[DV^]Va{Ku;)Vp{k\u;')/Z[0,0], 



where the [D'V] integration is over all velocity field con- 
figurations that are not related by Galilean invariance, 
and for each such field configuration the [i'Vc] integra- 
tion is over all field configurations related to the partic- 
ular field by a Galilean transformation. The expression 
Eq. ^ shows that using the ungauge-fixed functional 
Z, Eq. Il]), the calculation of a correlation function 
sums over configurations from all Galilean-related ref- 
erence frames, which in general will lead to an infinite, 
ill-defined quantity. This same situation would hold for 
a correlation function comprised of any combination of 
fields. Moreover, since these infinities are not just sim- 
ple terms that factor out, it means they can not be re- 
moved from the generating functional by simply dividing 
as Z[J,S]/Z[0,0]. 

The gauge- fixing term in Eq. ([5]) manifestly breaks the 
Galilean invariance of the action. We next demon- 
strate that this gauge fixed theory possesses a symme- 
try akin to the BRS symmetry of quantum field the- 
ory |12l . Il9l |20|, which when recognized, restores the 
Galilean invariance of the gauge fixed theory. We can 
multiply the generating functional Z'qp in Eq.([8]) by the 
constant / drjdr]* expjir;* • t]} where t] and rj* are con- 
stant complex conjugate Grassmann vectors {rji, rj*} — 0, 
T]^ = T]*^ = 0, and are independent of k and uj. For the 
"V^^-gauge" of Eq. ([5]) consider the following (infinitesi- 
mal) BRS transformation 



Sbrs'H = -7V0CC 



(^BRsVo = -cC(J7* 



V) 



SbrsV* 



-^VocC, 



(10) 

(11) 



where c is a constant with dimensions of velocity, ^ is 
a real Grassmann parameter (independent of k and uj) 
while all other modes of the velocity V(k, uj) and (T(k, uj) 
transform under (5b rs as implied by Eqs. ([3]) - (O, except 
we replace the velocity boost c by c ^ c({ri* + rj). The 
gauge-fixed action 

Sgf[V, t, rj, rj*] = S[V, a] + ^Vq ■Yo^iV*-V, (12) 

is BRS invariant: (5brS'5'gf[V, cr,r;,r;*] = 0. In the 
gauge fixed functional 



z 



GF 



(8) 



f [DY][Da]dridri* exp{-5Gi=^[V, a, r,, r/*] 
f dkduj{J ■ V + 1: ■ a) + e* -jy + e- 77*1(13) 



where and 9* are constant Grassman vector source 
terms, it is this BRS invariance that replaces the Gahlean 
transformation and leads to relations amongst correlation 
functions, which in analogy with quantum field theory, 
will be called the Slavnov- Taylor identities [l2, 120] ■ 

To derive these, start with the above gauge-fixed func- 
tional ZGi?[J, S, 0, 0*] and displace all fields by the in- 
finitesimal BRS transformation Eqs.([5][n]), with c -^ 
cC(*7* + ^7) a-nd Eqs. (|10llip . Since the measure and 
the gauge- fixed action Sgf are BRS-invariant, only the 
source terms will be displaced. Moreover, as (^ — 0, 
we can easily expand the exponential. The gauge-fixed 
functional transforms as Zqp — > Zgp + Sp rs Zqf , where 
^brsZgf = can be written as follows [18| : 



where 



Oj — j dkduj iJm{—)kj 
d S 



d S 

duj 6Jm{-) 



I ^rn\ jfSj 



dio (5S„(-) 



J,(-)<5(k)5H .(15) 



The shorthand (— ) stands for (— k, —u). 

We now show that the gauge fixed BRS symmetric 
theory Eq. P^ and the resulting Slavnov- Taylor identi- 
ties in Eq. (J14p are not simply a different way to write 
the known Ward identities 0] of this theory. First we 
demonstrate that the standard generating functional for 
the Navier Stokes equation Eqs. (H]) and ([2]), in the ab- 
sence of gauge fixing, contains spurious relations amongst 
correlation functions. To see this, return to Eqs. ()l|2p 
and perform the Galilean transformations Eqs. ^^■ 
Since the measure and action S are invariant under this 
particular change of variables, it leads to the relation 

(exp / dkdw {J(-) • V(k, uj) -f S(-) • (T(k, w)}) 

= (exp /dkdw{J(-) • V(k,w + c-k) 

- c-J(-),53(k),5(cc;)+S(~)-cr(k,cj + c-k)}), (16) 

where the averaging on both sides is with respect to the 
same measure and action S in Eq. ^. The above re- 
lation implies an infinite number of spurious relations. 
For example, acting with one derivative 5 / 5Ji{—\i.,uj) on 
both sides and setting J = S = 0, it gives (Vi(k, w)) — 
(Vi(k, w -I- c • k) — Ci(5'^(k)5(w)). Integrating over an in- 
finitesimal neighborhood near k = 0, w = 0, it leads to 
the relation c^ = 0, which is meaningless, since c is an ar- 
bitrary real vector which we are free to choose. Likewise, 
taking two derivative 6/5Ji{— 'ki, uji)5l6Jj{— 1^2,^2) of 
Eq. (J16p when integrated over an infinitesimal neigh- 
borhood near ki.2 = 0, loi^2 — leads to the relation 
-c^-AkiAwi(V,(b,0))-c,Ak2Aa;2(V^-(0,0))+c,Cj- =0, 
which again requires c = 0, which is senseless. In a 
similar fashion, any number of further derivatives of J 



and/or S fields will lead to an infinite number of spuri- 
ous relations. These examples show that at order c and 
all higher orders, the generating functional Eq. ([T|) con- 
tains an infinite number of spurious relations. Although 
many papers 0, d, [lO, [iJ, [3 have studied the generating 
functional Eq. ([T]), to our knowledge this Letter is the 
first to reveal this inherent inconsistency. 

At 0{c), an alternative way to reveal spurious relations 
is to work directly with the Ward identity expression. For 
the standard generating functional, the standard Ward 
identity relation found in the references [3, |^, llOj 13 
is OjZ — 0, where Oj is the operator in Eq. 



and Z is the generating functional Eq. ([T]). Consider 
(5/(5Ji(— k, —uj)Oj\3=-s=qZ = 0, this leads to the relation 
(fcj^Vi(k,cj) - %i53(k)(5(w)) = 0. Contract this with 
5ij, then fluid incompressibility implies that 6'^(\i)5{uj) = 
0, which is certainly false for k = and lj — Q. However, 
now we can exhibit an infinite number of such spurious 
relations all at 0{c). For example if two derivatives are 
taken (5/(5Ji(-ki, -wi)(5/(5Jm(-k, -Ci;)Oj|j=s=o2' = 0, 
then after contraction, fluid incompressibility implies the 
relation (5'^(ki)^(a;i)(Vm(k, w)) = 0, but the mean fluid 
velocity (V(k, tj)) can be generally nonzero for an ap- 
propriate noise force. So, this relation is also spurious. 
Moreover if one more derivative is taken with with re- 
spect to J or S, irrespective of the properties of the noise 
force, the resulting relation would be spurious. And sim- 
ilarly, any additional number of such derivatives will lead 
to an infinite number of spurious relations in which they 
will wrongly imply higher order velocity correlation func- 
tions are vanishing. 

In contrast, for the case of our gauge- fixed generat- 
ing functional Eq. (|13p . all spurious relations have been 
eliminated. At all orders higher than 0(c), this follows 
elegantly from the Grassmann property that C^ = 0. And 
at 0(c) the first nontrivial relation is obtained by differ- 
entiating, Eq. p4|) by S/SJoiS/S9j to give 



? 



{VoiVoj) + {vj*m) = 0, 



(17) 



which can be checked is correct with a similar relation 
found by replacing the differentiation by 9j with 9* . Fur- 
ther 0{c) relations can be obtained by taking more J 
and/or more 9 or 9* derivatives, but due to factoriza- 
tion, these are all proportional to the fundamental iden- 
tity Eq. nil), and so the GF-BRS functional Eq. ^ is 
completely free of spurious relations. 

These considerations support our statement that the 
generating functional Eq. ()13p is the correct and consis- 
tent one for this theory, and the correct expression of the 
Galilean invariance is the Slavnov- Taylor identities Eq. 
(fT4l) . found in this Letter. 

Up to now, we have been working with the generat- 
ing functional Zgf Eq. p3|) . but all our considerations 
can be applied to the generator of connected diagrams 
W = In Zqp and to the one-particle irreducible graphs 
r. Introducing the generating functional of one particle 



irreducible functions (i.e., the effective action) 

+ f dkdLj{J-V,i + 'E-(T,i), (18) 



we see that Jk — 



sr 



Sfc = 






dT 

Wf 



dT 

dnf 



The noise-averaged fields in the presence of the sources 
are denoted V^/ = (y) sources^ etc. In terms of T we can 
immediately write down the identity Eq. (|14p in the form 



= ^V^ 



6T 



^V^ 






ST 
drill 



(19) 



dv,i{\^,u) 



5T 



+ ki 



ST 



do; Sa-J,{-k,Lu) 



du) 



(5(k)(5(w) 



<5Kf(k,^) 
ST 



We now apply this formula to the problem at hand. The 
dependence of F on rj and rj* is simple since these con- 
stant Grassmann vector fields do not interact nor do 
they couple to the velocity or conjugate fields. Thus, 
we can immediately write r[Vc;, (Tc/, rye/, r/*;] = -i»7d ■ 
rici + r[Vc;,(Tc/], where r[Vc;,(Tc/] does not depend on 
either rid or ?7*;, with the first few terms written in [25J . 
The pertinent identity we seek is then obtained by 
inserting T into Eq. (fl9|) . and then differentiating this 
with respect to d/di-jf S/5V['\\<i,uo)S/Sa''^{~\<i,~uu) fol- 
lowed by setting rj'^' ~ r]*j = V^' = cr'^' = 0. Note that 
the terms in Eq. (jl9p depending on the gauge parameter 
£ do not contribute, and that we end up with identity 

ummm 






(1,1) 



(24), 



(-k,-a-;k,^) = r^;;^(O,0;-k 



-r2;k,w). 

(20) 

This well known identity, has been explicitly verified up 
to one-loop order [l], 0] ■ The GF-BRS theory preserves 
this identity as it should, but this theory is moreover 
free of all spurious relations. In addition, the origin of 



this identity is seen through the GF-BRS approach as 
arising in the limit of vanishing mean translational ve- 
locity, but completely independent of the properties of 
the fluctuating velocit y. T his conclusion, consistent with 
the assertions in [13l . Il4i] . is seen here in a systematic, 
mathematically consistent, derivation. 

The implications and consequences of Galilean invari- 
ance for Navier-Stokes fluid dynamics have been out- 
standing questions for some time now [i|, [j, 0, 0, H S [lO, 
VTVi [l3l ll4| • The formulation of the Navier Stokes equa- 
tion subject to random forcing as a classical stochastic 
field theory 0, 0, B lid, lli lij, 13 opens up the way 
to apply the systematic techniques originally developed 
for quantum fields [20] ■ This Letter has shown that a 
field theoretic approach provides clear understanding of 
the Galilean symmetry through its close relation to gauge 
fixing and BRS symmetry; the constraint singling out a 
unique reference frame leads to a global symmetry with 
anticommuting parameters first used in quantum gauge 
field theory [19|. 

In summary, we have discovered a new symmetry 
in Galilean invariant classical stochastic field theories, 
which provides a powerful new way to think about the 
analysis of these systems. In this Letter we focused on the 
Navier-Stokes equation but in broader terms our analy- 
sis applies also to other Galilean invariant theories such 
as the KPZ equation [26|, magnetohydrodynamics and 
the Burgers equation. An immediate and important con- 
sequence of our result is that, within the presently un- 
derstood generating functionals for these theories, there 
are spurious relations which our gauge fixing procedure 
corrects. Although the treatment in this Letter through- 
out assumed a stochastic noise force, formally our treat- 
ment is also applicable when the noise force is absent. 
Our work, apart from now providing the mathematically 
consistent expression for functionals of Galilean invariant 
stochastic field theories, would also be of importance in 
any numerical/simulation calculations of these theories. 

We thank W. D. McComb for useful discussions. Sup- 
port was provided to A.B. by UK PPARC and D.H. by 
CSIC-INTA (Spain). 
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